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Molecular dynamics study of a thin (one to five layers) lu- 
bricant film between two substrates in moving contact are 
performed using Langevin equations with an external damp- 
ing coefficient depending on distance and velocity of atoms 
relative the substrates, motivated by microscopic configura- 
tions. They show that the minimal friction coefficient is ob- 
tained for the solid-sliding regime. A detailed analysis of the 
results, the comparison with other microscopic modeling ap- 
proaches of friction, and the evaluation of quantities that can 
be compared to experiments, such as the velocity of the tran- 
sition from stick-slip to smooth sliding, are used to discuss 
the relevance of the microscopic simulations of friction. 



I. INTRODUCTION 

The problem of friction between two substrates which 
are in moving contact is very important technologically 
as well as very rich physically jL-f7). Following the de- 
velopment of atomic force microscopy, tribology has ap- 
proached the microscopic level and these studies are ex- 
panding because nanotechnology is now building devices 
that are so miniaturized that they begin to probe the 
microscopic properties of the materials. But the atomic 
friction microscope has its own difficulties, in particular 
for "high-speed friction" , because it does not operate at 
high velocities. This explains the interest of numerical 
simulations of friction || (see also the review paper |?]] 
and references therein). Moreover molecular dynamics 
can provide detailed information on the motion of indi- 
vidual atoms which help in the understanding of the basic 
mechanisms of friction, a problem which is still open in 
spite of the efforts devoted to it in the last few years. Un- 
fortunately, large-scale simulation of dynamical processes 
with three-dimensional models and realistic interatomic 
potentials are still too time consuming for large systems, 
mainly due to slow relaxation processes in such systems. 
Moreover the standard molecular dynamics methods con- 
sider the atomic degrees of freedom only and ignore the 
electronic ones, that may be a too crude approximation 
in modeling metal substrates. Therefore, one has to use 



simplified models such as two-dimensional f§-|ll| or even 
one-dimensional |l2|,[l3| ones. 

An intermediate complexity model, where a few 
atomic-layer film is placed between two rigid substrates, 
so that the top and bottom lubricant layers may stick 
to the substrates while leaving mobile atoms in between, 
was analyzed in a series of papers and also m 
our previous publications for commensurate [^p| and in- 
commensurate plf lubricant films. Thermal equilibrium 
was achieved with the help of Langevin equations with 
a random Gaussian force and a viscous damping with a 
constant coefficient rj. The Langevin thermostat was ap- 
plied either to all degrees of freedom of the system |^0|,|l| , 
or only to the ones perpendicular to the sliding direction 
p4] p"9[ . Owing to the relative simplicity of such models 
they allow the simulation of friction for a broad range 
of parameters and the study of general trends and laws 
of the phenomena. The main results obtained in these 
studies are the following. First, it was shown that the 
static frictional force f s for the contact of bare substrates 
is very small in most cases, except in the rare situation 
when the two surfaces are commensurate and perfectly 
aligned 0,^2|,^3|. However, even a sub-monolayer of mo- 
bile atoms between the surfaces leads to a finite f s which 
is proportional to the applied load and is almost indepen- 
dent on the system parameters [ fl4| . Second, a thin lubri- 
cant film (with a thickness smaller than ~ 10 molecular 
diameters) confined between two solids is always layered 
and often solidified because the confinement decreases 
the entropy of the film and shifts the bulk melting tran- 
sition to higher temperatures and lower pressures (see, 
e.g., and references therein). Third, when an exter- 
nal force / drives one of the substrates and it starts to 
move, the lubricant film melts and its width increases by 
~ 10% " 
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|l£j|. After the shear-induced first-order melt- 
ing transition, the lubricant exhibits a layer-over-layer 
sliding with strong two-dimensional order within shear- 
ing planes of atoms, where each layer moves coherently 
as a whole |T|,|2l]]. If then the driving force or velocity 
decreases back to a smaller value, the film solidifies again 
either in the solid state in the case of spherical molecules, 
or in an amorphous (glassy) phase for the case of lu- 
bricant consisting of long-chain (organic) molecules |L7j . 
Thus, if the top substrate is coupled through a spring to a 
stage moving with a constant velocity, then at the begin- 
ning the spring stretches and the force increases. Even- 
tually / exceeds / s , the top substrate begins to slide and 
the lubricant melts. Then the top substrate accelerates 
to catch up with the stage and / decreases, so the sub- 



1 



strates stick again. This periodic melting-freezing pro- 
cess was used by Thompson and Robbins [jTsf] for expla- 
nation of the stick-slip motion observed experimentally. 
Forth, when the film consists of flexible chain molecules, 
then it may be trapped in a glassy state, and the slid- 
ing occurs only at the lubricant / substrate interfaces (the 
plug- like flow) |i~7| . This explains why the effective vis- 
cosities of thin films may rise more than five orders of 
magnitude above bulk values (experimentally, the relax- 
ation times are fractions of a second for thin confined 
films and nanoseconds for a bulk). Finally, the regime of 
smooth sliding with atomic-scale velocities was studied in 
p8| , p0|j2l] . It was shown that the kinetic frictional force 
can be described by the introduction of a phenomeno- 
logical parameter which describes an "intrinsic" damping 
within the lubricant due to anharmonic coupling between 
its different modes. Also, He and Robbins in a recent pa- 
per JljJ studied the friction at low velocities. Their sim- 
ulations lead to a frictional force which is proportional 
to the load in accordance with the Amontons law, and 
rises logarithmically with velocity, again in agreement 
with experiments. The simulation results were explained 
as thermally activated motion of the lubricant over the 
periodic substrate potential. 

However, all these studies opened the question of the 
validity of results obtained with simple models, in partic- 
ular which results are of general validity and which ones 
are model dependent, and to what extend the results ob- 
tained on microscopic samples can be used to analyze ac- 
tual experiments performed on a macroscopic scale. We 
think that these questions are important for the future of 
molecular dynamics simulations of friction, and we want 
to address some aspects of them in the present work. 

A critical question is the energy flow out of the friction 
zone and the way it is described in the models because 
it governs the value of the friction coefficient. Thermal 
equilibrium is rather easy to simulate correctly because 
it does not depend on the value of the external damp- 
ing r\ in Langevin equations or on the method used to 
introduce the damping (although the rate of approach to 
equilibrium depends on r\). It is clear, however, that the 
non-equilibrium state of a system driven by an external 
dc force must depend on the method used to take into ac- 
count the energy exchange with the outside, for instance 
on the type of external damping which is assumed. One 
may predict that the dependence is not crucial if the rate 
of energy exchange between different modes of the sys- 
tem, which emerges due to nonlinearity of the equations 
of motion, is (much) larger than energy transfered outside 
of the system, described by the external damping 77. This 
case corresponds, for example, to a lubricant close to the 
melting temperature studied in |{li}-[l9|, especially for lu- 
bricants with complex (long-chain) molecules. But for a 
monatomic lubricant or a lubricant with simple molecules 
at low temperatures and low driving velocities, when the 
lubricant film is in the solid state, the simulation results 
may strongly depend on the external damping. In the 
present study we investigate this question by using a re- 



alistic energy loss mechanism. Comparing with models 
previously used in friction simulation (as discussed above, 
see also and references therein), the main new feature 
of the present study is that we use Langevin equations 
with an external damping coefficient n(z,v) which de- 
pends on the distance z of a lubricant atom from the 
substrate and on its velocity v with respect the substrate. 
Some significant differences with simpler modeling of the 
energy losses are exhibited. 

With this model, we pay a particular attention to the 
motion at small velocity «top of the top substrate, looking 
for the smallest possible value of the velocity for which 
the motion stays smooth. We show that the microscopic 
transition from smooth sliding to stick-slip motion takes 
place at a velocity which is many orders of magnitude 
higher than that observed in macroscopic experiments. 
This is an important point in the connection between 
simulations and experiments, which suggests that the 
macroscopic mechanism of the transition from stick-slip 
motion to smooth sliding is completely different from the 
microscopic one. 

We also used our extended model to examine various 
aspects of friction, which will be presented in the result 
and discussion sections. The paper is organized as fol- 
lows. First, in Sec. [0] we describe the model and the algo- 
rithm used in the simulations, and we introduce an exter- 
nal damping coefficient for Langevin equations, based on 
microscopic considerations, which attempts to describes 
the energy exchange between the lubricant atoms and the 
substrates in a realistic way. Then in Sec. [II we present 
the simulation results and analyze them considering the 
excitation of vibrations in the system. Finally, Sec. IV 
discusses the applicability of microscopic simulation to 
macroscopic experiments. 



II. MODEL 

We are interested in the regime of boundary lubrica- 
tion, when two sliding surfaces come in solid contact, 
and in the contact region, the surfaces are separated by 
at most a few layers of lubrication molecules. We con- 
sider a model with a thin (one to five layers) solid lubri- 
cant film between two solids. There are two reasons for 
this choice. First, there is always some lubricant (called 
"third bodies" by tribologists) between the surfaces. It 
may correspond either to a specially chosen lubricant, 
or to water or other material adsorbed from air, a dust, 
grease, wear debris produced by sliding, etc. In general 
the lubricant is not commensurate with the substrates, 
and, even if they are commensurate, the lubricant and 
substrates are seldom perfectly aligned. The external 
load squeezes the lubricant out from the contact area, so 
a very thin lubricant film is often formed. The last layer 
is however very difficult to remove |]J24|] so that there is 
at least one lubricant layer between the substrates. Due 
to the load, the lubricant is always strongly compressed 



2 



and thus has to form a closely packed structure. As a 
result, the confined thin film solidifies [QJitJ as was men- 
tioned above. Second, we are interested in low friction 
situations, and it is in the case of a solid thin lubricant 
film that we may expect to obtain the minimal friction. 
When the lubricant is liquid or amorphous, the kinetic 
friction is usually much larger p}-|6|, |is| , p5| . 

Moreover, contrary to our previous studies |pp| , pl| , now 
the surface layers of the substrates are simulated directly. 
The substrates consists of two parts. A rigid part forms 
the boundaries of the model system, and a deformable 
substrate layer is in contact with the lubricant, providing 
a much more realistic description of the substrate than 
previously. 
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FIG. 1. The model. 

Our three-dimensional system comprises a few atomic- 
layer film between two parallel rigid top and bottom sub- 
strate parts as shown in Fig. [j]. Each rigid substrate part 
has Ns = 132 atoms henceforth called s-atoms organized 
into a 12 x 11 lattice of square symmetry with the lattice 
constant a s = 3. The atoms of the bottom rigid substrate 
part are fixed while the top substrate part moves rigidly. 
Between the rigid substrate parts we insert atoms of two 
different kinds: 2N$ s-atoms model the surfaces of the 
substrate, which are therefore modeled more accurately 
than with rigid body, and N a iNi 1-atoms ("lubricant" 
atoms) form the lubricant film. In the x and y directions 
we use periodic boundary conditions. The top substrate 
part and all the atoms may move in all three dimensions. 
To shorten the expressions we henceforth call top and 
bottom "substrates" the rigid parts of the substrates. It 
should be kept in mind that these "substrates" are in fact 
only a part of the top and bottom solids that slide on each 
other. All atoms interact via a 6-12 Lennard-Jones (LJ) 
pairwise potential 



but the parameters of the potential (|l|) are different for 
different kinds of atoms. For the s-s interaction we took 
V ss = 3 and r ss = 3, for the 1-1 interaction, Vu = 1 
and Tu — 4.14 ("incommensurate" case), and for the 
s-1 interaction, Vi s — 1/3 (corresponding to a weak in- 
teraction between the substrate and the lubricant) and 
r is = \{ r ss +i"u) = 3.57, respectively. The cutting radius 
for the interaction was chosen as r* — 1.49 ru = 6.17. 
With these energy parameters, the s-atoms stick to the 
top and bottom substrates, covering them with mono- 
layers, while the 1-atoms fill the space between the sub- 
strates. This simulates a lubricant between two solids, 
and the surface (utmost) layer of each substrate is treated 
explicitly The equilibrium configuration of the lubricant 
corresponds to Ni layers, each having N a i = 80 atoms, 
organized into a close-packed triangular lattice slightly 
distorted by the substrate potentials. 

For the masses we took to; = m s = ms = 1 
which gives a typical frequency for the system of u s = 

1/2 

[V" s (r ss )/m s ] 1 = 4.9, and the corresponding char- 
acteristic period is r s = 2tt/lu s = 1.28. A typical 
frequency of vibrations within the lubricant is ui = 

1/2 

6 (2Vu/mirfA = 2.05. To each atom of the top rigid 
substrate we apply a dc force consisting of a driving force 
/ along the x axis and a loading force /i oat j = — 0.1 along 
the z direction. These parameters have been used in most 
of the simulations, except when otherwise specified. 

Units. Although we work with dimensionless quanti- 
ties, the numerical values of the model parameters have 
been chosen such that, if energy were measured in elec- 
tron volts and distances in Angstroms, we would have 
realistic values for a typical solid. However, for a dis- 
cussion of the applicability of simulations to real phys- 
ical systems, it is useful to couple the "natural units" 
(n.u.) used in the simulation with, e.g., the SI system of 
units. The basic parameters which were unchanged in all 
simulations, are the amplitude of interaction within the 
substrates (V ss = 3), which sets the energy parameter, 
the substrate lattice constant (a s = 3) which sets the 
length scale, and the mass of lubricant atoms (to/ = 1) 
as the mass parameter. A real system can be character- 
ized by the amplitude of interaction in the substrates V ss 
measured in eV, by the substrate lattice constant a s mea- 
sured in A, and by the mass of lubricant atoms fhi mea- 
sured in proton masses m p . If we introduce the following 
coefficients, v e = V ss /V ss = V ss /3, v r = a s /a s = 5 s /3, 
and v m — to//(100to/) = to;/100, then for a typical sys- 
tem we have v e ~ v r ~ v m ~ 1, and it is easy to find 
these coefficients for a real physical system. 

Now we can couple our natural units with the SI 
system of units. Namely, we have for the unit of 
length 1 to = lO 10 ^ 1 n.u., for the unit of mass 
1 kg = 6 • 10 24 n.u., for the unit of energy 1 J = 
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6.25 • lO 18 ^ 1 n.u., for the unit of force 1 N = 6.25 
10 8 (y r /vg) n.u., for the unit of pressure 1 N/ra 2 = 



6.25 ■ 10 (v^/Vf.) n.u., for the unit of time 1 s = 

l li 

0.98 • 10 13 \y & lv m iq\ n.u., and for the unit of ve- 
locity 1 m/s — 1.02 ■ 10~ 3 (^m/^e) 1 ^ 2 n.u. In particu- 
lar, the load force /i oa d = —0.1 n.u. used in the simu- 
lations, corresponds to the pressure P = — /load/a^ = 
1.11 • 10~ 2 n.u. = 1.78 • 10 9 N/m 2 . To compare with ex- 
perimentally used values, note that a realistic pressure is 
P ~ 10 7 N/m 2 , and the maximum pressure above which 
the plastic deformation begins, is P « 2 • 10 8 N/m 2 for 
gold (a minimal value for metals), P « 10 9 N/m 2 for 
steel, and P « 10 11 N/m 2 for diamond (the largest pos- 
sible value). As for velocities, a typical value when the 
transition from a stick-slip motion to smooth sliding is 
observed experimentally is v c ~ 1 fim/s — 10~ 9 n.u. [|4j. 

Equations of motion. We use Langevin equations for 
all mobile atoms, i.e. atoms that do not belong to a rigid 
substrate, 
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where ri a = {xi a . Via, zt a } is the coordinate of the zth 
atom and a — s or a — 1 for "substrate" or "lubricant" 
atoms respectively. 

The force /( mt ) is due to interaction between the mo- 
bile atoms in the system, 
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where the sum includes all "mobile" s- and 1-atoms except 
the given (ath) one. 

The last term in Eq. (Q) describes the interaction of a 
"mobile" s- or 1-atom with the bottom (S = 1) and top 
(S = 2) substrates. The force fi a s itself consists of three 
contributions as usual in Langevin equations, 



fia,S 



, (int) , (fric) /-(ran) 
Jia,S Jia,S "r J ia,S 



The first contribution f^g comes from the potential in- 
teraction of a given zth atom with all "immobile" atoms 
of the S'th (bottom or top) rigid substrate, 
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where the sum now includes all "immobile" s-atoms of 
the corresponding substrate and Ris is the coordinate of 
the ith atom of the 5*th rigid substrate. 

The second and third contributions in Eq. (Q) describe 
the energy exchange between mobile atoms and the rigid 
substrates. They approximately take into account the 



missing degrees of freedom of the substrates. More pre- 
cisely, the second contribution fj^^ describes a viscous 
damping when an atom moves relative the corresponding 
substrate, 
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where n(. . .) is the "external" damping coefficient which 
depends on the velocity and distance relative to the sub- 
strate, 

«(••■) = V (Zrel, Viei) , 

z IC \ = (-l) (s ~ lS, (z ia - Z S ), 1W = ha - Rs, (7) 

Rs = {As, Yg, Zs} is the center of mass coordinate of the 
5th substrate (for the bottom substrate we took R\ = 0). 
Finally, the third contribution Z^ 3 ^ in Eq. (^) describes 
the random (Gaussian) force acting on the ith atom from 
the S'th substrate. The amplitude of this force is deter- 
mined by the substrate temperature T, i.e. the corre- 
sponding correlation function is 
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The function rjR{z vc \) in Eq. (J 7 ^) coincides with the ex- 
ternal damping coefficient r](z Ie {) if the latter does not 
depend on the velocity (but may depend on the coordi- 
nate). Otherwise, if the external damping depends on 
f re i, the two coefficients are coupled by the relationship 



Vr(z, v, T) — / °° de e e n(z,v(e)), 



(9) 



For the top rigid substrate we use Newton equation of 
motion, 



M S R 2 = Nsf c: 



Fs, 
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where Mg = Ngmg is the mass of the rigid substrate, 
/cxt = {/, 0, /ioad} is the external force applied to it, 
and Fg = — fia,S=2 according to third Newton law 
(conservation of the total momentum of the system) . As 
we checked numerically, this technique leads to Gaussian 
distribution of velocities for all atoms as well as for the 
top rigid substrate with a correct width for a given tem- 
perature T. 

External friction coefficient. Fortunately, the motion 
of a single atom or a sub-monolayer film adsorbed on a 
crystal surface, namely their vibration near the equilib- 
rium position, has been well studied experimentally and 
theoretically J2^j2^] . This allows us to model the exter- 
nal damping of lubricant atoms near the substrates with 
a reasonable accuracy. For the external damping coef- 
ficient r\(z, v) which models the energy loss of an atom 
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into the substrates and enters into Eqs. (|6[g), we take 
into account its dependence on the distance z from the 
corresponding substrate and on the relative velocity v 
according to r](z,v) — r]i(z)r]2(v). First, we assume that 
the damping rate exponentially decreases when an atom 
moves away from the substrate, and saturates at some 
level when the atom approaches very close to the sub- 
strate. 



771(2) = 1 — tanh[(z — z*)/z* 



(11) 



The characteristic distance z* is chosen as the height of 
the pyramid with the base constructed of four rigid sub- 
strate atoms with square length a s = 3, and the vertex 
with the s-atom at the distance r ss = 3 from the rigid- 
substrate atoms; this leads to z* = 2.12. Thus, for the 
atoms in the s-layer, where z ~ z*, we have 771 ~ 1, 
while for the atoms in the first (closest to the substrate) 
lubricant layer we obtain 771 ~ 0.1. 

Second, to determine 772 (v) we used the results known 
for the damping of the vibrations of an atom adsorbed on 
the crystal surface. According to the theory [^9 30 1 , when 

ibrium 



an atom vibrates with a frequency ui near its equi 
position, the oscillation decays due to creation of phonons 
in the substrate with the rate 
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where the surface local density of phonon states may be 
described approximately by the function [B0| 
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and uj m is the maximum (Debye) phonon frequency of 
the solid substrate. The one-phonon damping mechanism 
works for frequencies lu < u> m . At higher frequencies the 
damping is due to multi-phonon processes. Moreover, in 
the case of metal or semiconductor substrates there is an 
additional damping due to the creation of electron-hole 
pairs in the substrate so that the corresponding damping 
coefficient is typically of order 77 ~ 10 _2 o; m (see p9| , [30| ]). 
Using these results, we took for the dependence r]2(v) the 
following approximate expression, 



m(v) = Vmin + T] ph (2nv/a), 



(14) 



with a = a s for the motion along the substrate and a = z* 
for the motion in the z direction. For the minimal contri- 
bution ?7 m i n we chose ?7 m in = 0.01 uj s — 0.049. This value 
is compatible with the justifications given above, and it 
leads to a thermalization of the system in reasonable sim- 
ulation times t <3 -10 3 r s . The cutoff (Debye) frequency 
was found from the phonon spectrum (the Fourier trans- 
form of velocities of all atoms and the top rigid substrate 
calculated for / = at T — 0.025) which gives oj m = 15. 
Thus, as a function of frequency the external damping 
behaves as 
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so that at small frequencies, to <C to m , we have 772 (to) ~ 
Vmin + (16A4J w 4 « (4.9 + 0.47w 4 )-10~ 2 , while at the 
frequency u> = (4/7) 1 / 2 w m = 11.34 the external damping 
reaches its maximum value 772 ~ 22. 

Energy losses. When one part of a system moves with 
respect to another part with a relative velocity v, then 
the rate at which work is done is equal to e — vf, where 
/ is the total force acting on the former. This may be 
used to define the energy losses of a given atom as 
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where the sum is over all atoms of the system including 
the s-atoms of the rigid substrates, and fw is the force 
acting on the ith atoms from the i'th one. Note that this 
force has to include the damping and random (Gaussian) 
contributions of the total force when i' corresponds to the 
s-atom of rigid substrates. Then, taking a corresponding 
sum over atoms and averaging over time, one can find the 
energy losses in a given atomic layer of the lubricant, or 
those in the substrates, as well as separate contributions 
from different degrees of freedom. 

To illustrate this definition, let us consider an ideal 
case, when the bottom substrate is immobile, the top 
substrate moves with a velocity v in the x direction, 
and the lubricant film with the velocity while the 
force F is applied to the top substrate. We have 

= -Fbot/lub = — -Flub/bot = F, 



-Flub/top — — -Ftop/lub 



and we obtain e top = — i(— F) (v — |w) = jFv, £i u b = 

-\ \ F (!« - «) - F M = £ bot = -\f (0 - \v) = 

jFv, thus the total losses are equal to e = Fv which 
coincides with the work done by the external force as 
expected. 

Curved substrates. In most simulations both top and 
bottom substrates were flat. However, we have done also 
a few runs with curved substrates ]3l| ] , where the z coor- 
dinate of the bottom rigid substrate is defined by 



(17) 



where L x y is the size of the system in the x or y direc- 
tion and h x ^ y are the corresponding curvature parame- 
ters. Similarly the z coordinate of the top rigid substrate 
is defined by the expression 
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where { X2 , Y2 , Z2 } are the center of mass coordinates 
of the top rigid substrate. 

Algorithm. In most of the simulations we used the 
constant-force algorithm [^o|,^l). Namely, at the begin- 
ning we equilibrate the system at a temperature T (in 
this paper we present the results for T = and T — 0.025 
which corresponds to room temperature for energies mea- 
sured in eV). Then we adiabatically increase the shear 
force / in steps A/ (typically A/ = 10~ 3 or smaller). 
Each increase of / is accomplished progressively in i t 
(typically i t = 200 or more) substeps each of duration 
t s . When / has reached its new value, we wait dur- 
ing a time i t r s in order to allow the system to reach a 
stationary state, and then measure the atomic coordi- 
nates and velocities during the next time interval itT a . 
Thus, the force is changing with the average rate 1Z/3t s , 
where 1Z — Af/it. In some simulations we also used 
the constant-velocity algorithm (which is typically used 
in tribology simulations, see &), when a constraint is 
used to keep constant the top substrate velocity. We 
found, however, that the const-w algorithm leads to a 
much smaller accuracy of the results than the const-/ 
one, probably because the friction force in this case is 
not a self-averaging value. Thus, here we present the 
results of the const-/ runs only. 

Most of the simulations have been done for a rather 
small system. The issue of system size raises the usual 
question in computer simulations of whether it is bet- 
ter to take a small system and study it carefully or to 
take a larger system and do a less comprehensive study 
over shorter times. Our simulations show that the pro- 
cesses under investigation are characterized by long re- 
laxation times. The approach with adiabatic changing 
of the driving force necessitates going through a whole 
cycle of / changes. For a much larger system it is not 
possible to reach equilibrium in a reasonable time for 
this type of study. However, we also made few runs 
for a system four times larger and found that the re- 
sults did not change. Besides, in physical systems mi- 
croscopic contacts are of a size similar to that of our 
system, for example 10 x 10 atoms. Even in specially 
prepared surfaces, with experiments performed in high 
vacuum conditions, ideal surfaces usually do not extend 
beyond about 30 lattice constants. Of course the results 
deduced from such simulations of ideal surfaces should 
not be used without extreme caution to analyze the re- 
sults of macroscopic experiments. We shall come back to 
this point in the discussion because the connection with 
macroscopic experiments is a difficult point for all stud- 
ies of microscopic friction, whether they are experimental 
(with atomic force microscope) or numerical. 



III. RESULTS 

Before going into the details of simulation results, let 
us discuss some general features. Typical dependences 
of the velocity of the top substrate in the x direction, 
^top(/)i and the change of its z coordinate, Az top (/), as 
functions of the driving force / for the one- and five-layer 
lubricant films are shown in Fig. |2| Each plotted point is 
the average over i t = 200 data points recorded at times 
separated by At — r s and corresponds to the steady-state 
motion for a given /. 




10" 3 10" 2 10"' 10' 
force 




10" 3 10" 2 10"' 10° 

force 

FIG. 2. (a) The velocity of the top substrate and (b) 
the change of the lubricant thickness for the one-layer (solid 
curves and open diamonds) and five-layer (dotted curves and 
solid diamonds) lubricant films. 

Let us analyze first the behavior of the five-layer film. 
The motion starts when the force exceeds a static friction 
force f s sa 0.0017—0.0018. One can see a plateau on the 
v top(f) dependence for the force interval 2-10~ 3 < / < 
0.03, then at larger forces v iop approximately linearly in- 
creases with /, and finally at / > // w 0.30 — 0.31 the 
"regular" motion is destroyed. The velocity of the top 
substrate for these force intervals lies within the interval 
1 < Vtop < 4.5, so that the washboard frequency, defined 
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as w was h = 7r(w t op)/a, lies just inside the phonon spec- 
trum of the lubricant film as will be discussed in details 
later. In this force range, the width z of the film increases 
as shown on Fig. ||b. The system exhibits hysteresis: if 
one starts with a force / < // and then decreases it adi- 
abatically slowly, the system comes back to the T = 
ground state for a force /& ~ 2 • 10~ 5 , which is much 
lower than f s . The velocity of the top substrate at the 
backward transition drops down abruptly from a value 
Vb w 0.03 to zero. Moreover, if we start from the sliding 
state and then remove the driving force at all, the system 
again comes to the T = ground state. 

The behavior of thinner lubricant films is very similar. 
For example, the one-layer film starts to move at f s sa 
0.019 — 0.02. The forward transition now takes place at 
// w 0.43 which is slightly larger than that for the five- 
layer film, although the velocity at which the film starts 
to melt is approximately the same, vj ~ 4. The backward 
transition takes place at w 1.5- 10 -4 when vi w 0.085. 




FIG. 3. The static friction force f a as function of the 
number of lubricant layers TV; for T = 0.025. 

Although the sliding regime is approximately the same 
for different lubricant widths, the static frictional force 
needed to initiate the sliding decreases approximately ex- 
ponentially as the width increases as shown in Fig. ||. The 
dependence f s (Ni) may be explained qualitatively in the 
following way. To initiate the sliding the lubricant atoms 
must overcome barriers created by the substrate atoms. 
But the height of the barriers strongly (exponentially) 
depends on the distance z of the lubricant atoms from 
the substrate. It decreases when the lubricant atoms are 
moved away from the substrate. In the one-layer case 
the confined lubricant has to adjust to both substrates 
simultaneously and thus has no freedom to move in the 



z direction. On the other hand, when the lubricant is 
thicker, its utmost layers can move in the z direction due 
to elasticity of the film, and thus can allow the atomic 
positions to find a minimum-energy saddle configuration. 
This is confirmed by the analysis of the system configu- 
rations just before f s . They show that the shifts of the 
atoms in the z direction from their / = equilibrium po- 
sitions in the five-layer film are much larger than those 
in the three-layer film which in turn are larger than in 
the two-layer film. 

The dependence of the static force f s on the tempera- 
ture T is only significant for Ni > 1 when the lubricant 
atoms have enough freedom to move in the z direction. 
For example, for the five-layer film we found f s « 0.0017 
at T = but f s m 0.001 at T = 0.025. 




(b) 
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1 2 3 4 5 

10" 3 time 

FIG. 4. Beginning of the sliding of the five-layer film, 
(a) The x- velocity of the top substrate (solid curve) and the 
lubricant (dashed curve). Inset: the same enlarged for shorter 
times, (b) z t o P versus time. 

The initial stage of the sliding is illustrated on Fig. 
||, where we plot v x (t) for the top substrate and the 
lubricant and ztop(t) for the five-layer film. One can 
see that at the beginning, the top substrate (the rigid 
substrate with one attached s-layer) starts to move 
first. Its velocity increases according to a law v tov {t) ~ 
[f s /(m 3 +m S )Ves{t)] [l-e-wOHl], where the ef- 
fective friction coefficient is defined as ^ c ff = f/( m s + 
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ms)vt op , and achieves a plateau in the steady state. Soon 
after, however, the lubricant is also involved into the 
motion and reaches the velocity ui UD = |«top(/s) m the 
steady state. Simultaneously the width of the lubricant 
slightly increases. We observed that for our system sizes, 
the lubricant begins to move as a whole. All its atoms 
begin to move almost simultaneously. However, for larger 
sizes of the contact one could expect that the sliding may 
begin with the creation of "moving islands" [ ^2[ . 

An interesting result is that the one-layer film can be 
in the sliding steady state for much larger dc forces than 
thicker films. The solid sliding exists for forces within the 
interval 1.1 < / < 3.7 when the velocity is 5 < v top < 9. 
This steady state cannot be reached with by adiabatically 
increasing the force because the lubricant melts earlier. 
But it can be obtained with sharp increase of the force, 
e.g. if one takes the / = 0.2 steady state and applies the 
force / = 2. This high-speed steady state corresponds to 
the "flying" regime predicted in j2l]] and will be described 
in more detail below in Sec. [V. 



A. Solid sliding 

All steady states described above correspond to a solid 
sliding regime, in which the lubricant moves as a whole 
with a velocity equal to half of the velocity of top sub- 
strate. The distributions of velocities for all forces can be 
well approximated by Gaussian curves if we use different 
"temperatures" for the lubricant and the s-atomic sub- 
strate lavers as well as for different degrees of freedom. 




force 

FIG. 5. Lubricant temperatures versus force for a 
five-layer film. Inset: velocity histograms (dots) and the cor- 
responding Gaussian distributions (solid curves) for / = 0.02. 

For example, for the five-layer film we show in the in- 
set of Fig. u the velocity histograms for the / = 0.02 



case by dotted curves and the corresponding Gaussians 
by solid curves. Thus, we can define the "temperature" 
for a given degree of freedom (a = x, y, z) as well as 
for a given layer I of the lubricant or the substrate as 
T a = m((v a — (v a )) 2 ), where (. . .) designates the aver- 
aging over time and over all atoms in the given layer, and 
then use the values T a (l) (let / = corresponds to the s- 
atomic layers of the bottom substrate, I = 1 corresponds 
to the first layer of the lubricant film, etc.) as reliable 
characteristics of the driven system. The distribution of 
temDeratures across the svstem is shown in Fis. 




-1 
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FIG. 



Distribution of (a) the temperature T z and (b) 
the energy losses across the five-layer lubricant for different 
dc forces / at T su b = 0. The numbers -1 and 7 correspond to 
the rigid substrates, the numbers and 6 correspond to the 
s-atomic substrate layers, and the numbers 1 to 5 correspond 
to the lubricant layers. 

First, note that comparing with the lubricant tem- 
peratures, the temperatures of s-layers of the substrates 
are very small (for the T = simulation). Second, at 
low forces/ velocities the temperature is not uniformly 
distributed over the lubricant film, the boundary layers 
which are in moving contact with the substrates, have 
a higher temperature than those in a middle of the lu- 
bricant. But at large forces/velocities, / ~ 0.02 — 0.2 
when i>top ~ 1 — 4, the lubricant temperature is approxi- 
mately uniform across the lubricant. This indicates that 
anharmonicity effects, which are responsible for energy 
exchange between different layers within the lubricant, 
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become large enough at these forces. Third, the lubri- 
cant temperature increases with growing / (see Fig. |^) 
until it finally melts at / = //. Forth, for all studied 
cases we found that T z ^> T x > T y . This indicates that 
the driven system is strongly out of equilibrium. 

The calculation of energy losses E\ oss for different 
dc forces shows that total losses are close to the ex- 
pected values. For example, for / = 0.002 we obtained 
E loss « 0.24, while fvtop « 0.26; for / = 0.02, E loss a 3.2 
(fvto P « 3.6); and for / = 0.2, E loss a 102 (fv top « 103) 
correspondingly. The energy is lost mainly due to the 
motion of atoms along the direction x of the driving, the 
y and z components of the losses are negligible. The 
distribution of energy losses in the normal direction of 
the system is shown in Fig. ||b. One can see that the 
energy is lost mainly within the rigid substrates and in 
the utmost lubricant layers (i.e. in the layers which are 
in moving contact with the substrates). 



F a (u) = 



E 

3 



dtvj a (t)e % 



(19) 




where the sum may be over all atoms or only over a se- 
lected subset, such as the lubricant atoms. The spectra 
calculated in this way for different degrees of freedom 
(a = x,y, z) separately for the substrates and the lubri- 
cant are shown in Fig. fj] for Ni — 1 and 5. One can 
see that the lubricant spectrum occupies the frequencies 
1 < uj < 6, while the substrate modes lie mainly within 
the interval 6 < uj < uj rn . The x and y spectra are broader 
than the z one. For the one-layer lubricant, the spectrum 
F z {ui) for the lubricant has a maximum at uj ~ 3 which 
is quite sharp. 

Calculating in the same manner the spectrum in the 
solid-sliding regime, we observed the excitation of z- 
oscillations of the lubricant with the frequency uj = cj was h 
as indicated by arrows in Fig. 0. We also found that (a) 
the intrinsic energy losses in the lubricant become im- 
portant only when w W ash is higher than the lower bound- 
ary of the lubricant phonon spectrum, w was h > 1; (b) at 
/ = 2-10 -3 one can see many higher-frequency harmonics 
of w W ash, i.e. the vibration of lubricant is highly anhar- 
monic; when the force increases even more, the energy 
exchange between the modes of the lubricant distributes 
the energy among all the lubricant modes; (c) at / = 0.02 
and / = 0.2 the lubricant is strongly heated, the energy 
of the translational motion is distributed over all three 
(x, y and z) degrees of freedom. 
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FIG. 7. Phonon spectra of the lubricant (solid curves) and 
the substrates (dotted curves) for one- and five-layer films. 
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To analyze in details the energy losses for the solid- 
sliding regime, we have to know the phonon spectrum of 
the system. To find it, we saved the coordinates of all 
atoms at / = and temperature T = 0.025, and then 
calculated the Fourier transform 



FIG. 8. Change of the z-spectrum (|l9p of five-layer lubri- 
cant with increasing driving force. The arrows indicate the 
position of the washboard frequency. The curves for / 7^ 
are artificially shifted upstair. 
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The behavior of the thin Ni = 1 film is similar. 
The temperatures of the lubricant begin to increase at 
/ > 410~ 3 when the velocity is v top > 1 so that the wash- 
board frequency penetrates into the lubricant phonon 
spectrum. Then at / > 0.1 the excitation of lubricant 
vibrations strongly increases, especially for the z coordi- 
nate, so that the width of the film begins to increase (see 
Fig. ^d). The temperature reaches the value T z ~ 0.4, 
and soon the lubricant melts. On the other hand, for the 
high-force sliding state (the "flying" regime) the pertur- 
bation of the film is much smaller, Az < 0.5 and T z ^0.1. 
At small forces/ velocities the energy losses are mainly in 
the rigid substrates. At / ~ 0.02 — 0.2 the losses are 
approximately equally distributed between the two rigid 
substrates and the lubricant. Finally, in the high- force 
"flying" regime the main losses are in the lubricant film. 



B. Effective friction coefficient 

In a steady-sliding regime, where the top rigid sub- 
strate and one s-layer attached to it slide over the lu- 
bricant film, it is convenient to introduce the effective 
friction coefficient as 



T?eff(utop) 



(ms + m s )v iop (f) 



(20) 



Note that a change of the velocity of the top substrate 
can be expected to take a time of the order r rc i a x = %S ■ 
Using the results of simulations, one can see that at small 
forces the effective friction can be smaller than r] e s < 
10 -3 , so the relaxation times are larger than T re i a x > 10 3 . 

Energy loss can only occur through the externally im- 
posed damping. One part comes from the motion of the 
lubricant as a whole. This is the losses that one would 
get if the lubricant were sliding as a perfect solid rela- 
tively the rigid substrates. The corresponding contribu- 
tion to the effective friction coefficient is due to the exter- 
nal damping [with the coefficient r)(z, v) introduced into 
the equations of motion from the beginning] of the uni- 
form a;-motion of lubricant atoms relative the substrates 
with the velocity g^topj and yields a "universal" depen- 
dence ?4op (/) wn i cn does not depend neither on the 
lubricant width nor on the masses m s and ms- In what 
follows we call this contribution as the "perfect-sliding" 
one. A second part comes from the internal motions in- 
side the lubricant and the surface layers of the substrates. 
The driving excites atomic vibrations (mainly in the z di- 
rection as mentioned above) which then are distributed 
over other (x and y) degrees of freedom and finally are 
damped again due to the external damping i](z, v). This 
contribution will be called below as the "internal-losses" 



In the solid-sliding regime, when the top rigid substrate 
with one attached s-layer moves as a whole with a velocity 



( «top) j the bottom rigid substrate with one attached s- 
layer does not move at all, and the lubricant film moves as 



a whole with the velocity v\ u b 
frequency is equal to 



{Vtop) 



the washboard 



^wash = 27rw lub /a = n{v top )/a, 



(21) 



where one should take a = a s = 3. Let us write the 
balance of forces for the top substrate with one attached 
s-layer as 



F = Nsf = N aimir)* (v top ) V\ ah , 



(22) 



where we have introduced the total viscous damping coef- 
ficient ij*(vtop) for an atom in the utmost lubricant layer. 
Using the definition ( pp| ) of the effective friction coeffi- 
cient, we obtain the following relationship between the 
coefficients T] cS (v t0 p) and ry*(« to p), 



7?eff(wtop) 



mi 



N, 



a I 



l 



2 ms + m s N s 



V*(vtop)- 



(23) 



In a general case the total damping r]*(vtop) defined by 
Eq. (|22|), has to consist of two contributions, the "perfect- 
sliding" contribution ry e xt(^top) and the intrinsic losses 
f?int(ftop)- In the two following subsection we evaluate 
successively these two components of the losses. 



1. "Perfect-sliding" contribution to the energy loss 

In the perfect-sliding approximation the atoms in the 
utmost lubricant layers feel only the external damping 
??cxt(vtop) ~ m( z iubr) ?72(w W ash) due to energy exchange 
with the rigid substrates. Substituting the washboard 
frequency ( pl| ) into Eq. (^), we obtain 



7?ext(v t op) 



1 — tanh 



167^ 



J top 



1 - 



"top 



3/2* 



(24) 



Assuming that all the damping within the lubricant is 
due to the external one, r]*(vto P ) = rj ext (vt p) , we obtain 
the "universal" ("perfect-sliding") dependence 



(/) = 



2N, 



f 



N a i mir] eKt (vtop) ' 



(25) 



which does not depend neither on the number of lubricant 
layers not on the substrate masses m s and ms, because it 
corresponds to the steady state (while a delay of response 
of D t0 p when / varies non-adiabatically, has to depend on 
the masses). 

In particular, for the parameters used in simulation 
(mi — m s — ms — 1, Ns — 132 and N a i = 80), we obtain 
rj cS — 0.15 77*. Then, using in Eq. fl24| ) the values z* = 
2.12, tu m = 15, a s — 3, and taking from the simulation 
ziubr ~ 5.25 for the T = ground state configuration, 
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we obtain rj cx t(v) 



4.9 + 0.57 v 4 (l-0.49-l(TV) 3/2 

1CP 3 , so that the external phonon damping exceeds the 
minimal one (which models the electron-hole damping) at 
v top ^1-7 and reaches its maximum value 77 cx t ~ 2.16 at 
Vtop ~ 10.8. If all the damping within the lubricant were 
due to the external one, the effective friction coefficient 
would be equal to 



^ff ni) («to P ) - 0.728 -10- 3 [1- 



0.H6 vt op (1 



0.49- 10" 2 - 2 



,3/2 



top / 



(26) 



with the maximum r] c ff w 0.32. The "perfect-sliding" 
dependence ( p5| ) is presented in Fig. [| One can see that 
it agrees rather well with the simulation data at small 
(/ < 10~ 3 ) as well as at high (/ > 1) forces, when the 
washboard frequency lies outside the lubricant phonon 
spectrum and, thus, the internal motions of the lubricant 
are not excited. 
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FIG. 9. The "perfect-sliding" dependence ( |25| ) and the 
T = simulation results. 



An important result of this analysis is that the lubri- 
cant cannot support a force larger than / max ~ 7 when 
^top ~ 10.8. At larger forces/ velocities the system can- 
not dissipate the energy injected into it due to driving, 
and the solid-sliding regime must be destroyed. 

However, in the simulation we found that t>top(/) < 

v tap\f) an( l // < /max in all cases. To describe simu- 
lation results more accurately, the following three fac- 
tors have to be taken into account: (i) at velocities 
vtop > 1 — 2 the lubricant width increases (due to anhar- 
monicity of the LJ potential), (ii) there are additional 
intrinsic energy losses within the lubricant as well as in 
the s-layers, and (Hi) due to these intrinsic energy losses, 



the lubricant is heated during the sliding. The first factor 
increase of v top while the second and third ones tend on 
the contrary to decrease it so that they must be studied 
in details before a conclusion can be drawn. 



2. Excitation of lubricant vibrations 

When the lubricant film moves over the bottom sub- 
strate with the velocity i>i u b dragged by the top substrate 
with the dc force F = N$f, the dragging force per lu- 
bricant atom in the utmost lubricant layer is equal to 
/lub = (Ns/N a i)f. The total energy injected into the 
lubricant by the external dc force per time unit per one 
lubricant atom, is equal to 



£tot = /lub^lub = m lV* v lubi 



(27) 



where we used the definition (|22j ) of the coefficient r\* . 
This energy must be absorbed in the system. One chan- 
nel of energy dissipation, the energy exchange with the 
rigid substrates, was already considered in the previous 
subsection, and it yields 



£ext = TO/77 cx t«i ub . 



(28) 



An extra energy losses, £j nt = e to t — £cxt, must be at- 
tributed to intrinsic losses within the system. 

To study the intrinsic losses, we will use the technique 
similar to that described in |3^Jl^,^]. Let us consider 
the internal vibrations of the system as a damped oscil- 
lator £(t) of mass m, internal frequency too and damping 
coefficient 7y , where 770 corresponds to the so-called full 
width at half maximum (FWHM) of the spectrum (a 
generalization to many oscillators is trivial). If the oscil- 
lator is excited by an external force oscillating with the 
frequency u> and the amplitude /o, 



rnC + m m £ + rru4i = /(*) = Re (jbe^) , 



(29) 



the steady-state motion of the oscillator corresponds to 
its vibration with the frequency of the external oscillating 
force, 



£(t) = k Re 
m 



tLorio 



(30) 



The rate of energy losses (the energy absorbed by the 
oscillator per one time unit) is equal to 



w ,7?o,/o) = - 
r 



2 m 



(31) 



where r = 2ir/u> and 



p(u>) = Im 



1 



Ljri 



2\2 



(32) 
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In the solid-sliding regime the oscillating force is due 
to motion of the lubricant film over the corrugated sub- 
strates and thus is characterized by the frequency u> = 
<^wash- Then, assuming that the system can be described 
as k different oscillators and using Eqs. (p7| f32|) , we fi- 
nally obtain i]*(v top ) = r) CKt (v t o P ) + ?7int(wto P ), where 

4 k 

fftnt(Wtop) = — V^osc^wash; W$ , T)^ , ffi ) . (33) 

m l v top j=1 
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77o %>h(k>)/?7ph f J i wri ere jJq is a constant. Then, 
from the spectra of Fig. [j] we choose the following five 
frequencies: Uq — 7.3 which can be associated with the 

(2) 

z vibrations of the substrates, uiq = 5 which describes 

(3) 

the lowest x and y modes of the substrates, u = 3.3 
which corresponds to the lubricant oscillations, the fre- 
quency uJq 4 ^ = which describes the excitation of 
the lubricant by the second harmonic of w was h, and fi- 
nally (for the Ni = 5 case only) ui^ — 0.4 which may 
be associated with the vibration of the top substrate as 
a whole. The results of fitting are presented in Fig. nu, 
where the values of the fitting parameters 7^*' and ftp 
are given in the caption to the figure. One can see that 
this simple approach gives a reasonable agreement with 
the simulation data and shows that the velocity of the top 
substrate is determined by the "perfect-sliding" approx- 
imation at forces / < 1CP 3 , by the intrinsic losses cor- 
responded to excitation of lubricant vibrations at forces 
1CP 3 < / < 0.1, and by excitation of the substrate vi- 
brations at higher forces. However, we have to note that 
a general theory of kinetic friction which could give, e.g., 

the values of the parameters rj^ and /q from first prin- 
ciples, is still lacking. 
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FIG. 10. Fitting of the simulation results for (a) Ni — 1 
and (b) iV; = 5. Open diamonds are for the simulation data, 
the dashed curve shows the "universal" dependence (p5|), the 
dotted curve is obtained with Eqs. (|31r|33|), an d the solid 
curve is obtained from the dotted one by removing the unsta- 
ble branches where dvt op {,f) /df < 0. For the fitting we used 
the following frequencies, u^P = 7.3, lo^ = 5, lo^ = 3.3, 
uj^p = 1.65, Wq = 0.4, and the following fitting parame- 
ters for the widths of the resonance modes and the oscillating 
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(i) 



(4) 
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To fit the simulation data, we also assumed that 
the oscillator's damping coefficient depends on the 
frequency according to the expression tJq 1 (lu) = 



C. Minimal sliding velocity 

As mentioned above, the system exhibits hysteresis: 
when the force decreases starting from the sliding state, 
a backward transition takes place at / = /(,< / s , and 
the velocity drops down from a finite value v% op = vi, to 
zero. The same drop of the velocity is observed in free 
runs, when one starts with the solid-sliding steady state 
and then remove the driving: the velocity ft op slowly 
decreases till it reaches the minimal value Vb and then 
drops to zero as shown in Fig. O. It is well known (e.g., 
see |54|) that a similar hysteresis is observed for a sin- 
gle driven particle at T = 0. If a particle of mass M 
is placed into a sinusoidal external potential of (total) 
height £ = maxV(x) — minima;) and period a, and is 
driven by the dc force F, the forward locked-to-running 
transition takes place at F — F s = ir£/a, and the back- 
ward transition, at F = F b = (2\/2/ V) tjVMS, when 
the velocity is ~ y/aF s /M. In the underdamped 
case such that r\ < r\ c = {it 1 / '2v2j \/£/Ma 2 , we have 
Ft < F s . Therefore the system exhibits hysteresis due to 
the inertia of the particle. 

If the top substrate with the attached s-layer would be 
treated as a rigid object, taking a — 3 and mj + m s = 2, 
the critical damping would be 



fs 



a(m s 



1/2 



0.8 



(34) 



The effective friction coefficient just above the transition 
to sliding is r] 8 — f s /(ms + m s )v s . Hysteresis is expected 
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if T] s < j] c , i.e. v s > v bR , where 



v bR = -a 



7r 3 mj + m t 



1/2 



0.62 V/, 



(35) 



In the simulations the transition to sliding occurs for f a ~ 
10~ 3 — 10~ 2 while v s ~ 1, so, for the chosen set of param- 
eters, the system under investigation is clearly in the un- 
derdamped limit. We expect the backward transition at 
f b = ^v^AVtt) VcS \/( m s + m s ) af s ps 1.2477 e ffv77- 




time 

FIG. 11. The x velocities of the top substrate and all 
layers (including the s-layers) versus time for the "free" run, 
when the five-layer system starts from the steady state cor- 
responded to / = 0.2 and then the driving force is removed. 
Inset: enlarged part of the dependences at late times. 

Now we can check whether the condition v bR ps 
0.62 \fj~s agrees with the simulation data. For the solid- 
sliding regime with Vu = 1 at T = we found from free 
runs that for Ni = 1, f s ps (1.9 - 2)-10~ 2 , v b < 8.5-10~ 2 , 
while v bR ps 8.7-10- 2 ; for N t = 2, f s ps (8.9 - 9.1) 
< 6 • 10~ 2 , while v bR ps 5.9 • 10~ 2 ; for N t 
■10~ 3 , v b < 5-10~ 2 , while v bR ps 4 
f s ps (1.7-L8)-10 -3 , v b < 3-10- 2 , 
We also made runs for a smaller value of 
the interaction parameter Vu = 1/3, when the lubricant 
is weaker, and found that for Ni = 1, f s ps (1.3— 1.4)10 -2 , 
v b < 8 • 10~ 2 , while w 6fl ps 7.2 • 10" 2 ; and for A, = 5, 
/ s ps (1 - 1.1)-10- 3 , v b ~ 1.5-10- 2 , while v bR Ps 2-10~ 2 . 

Thus, this simple inertia consideration quite well de- 
scribes the hysteresis of the solid lubricant film and, at 
the same time, it shows that the minimal velocity of 
the top substrate in the solid-sliding regime cannot be 
smaller than v b ~ 10 -2 n.u. ~ 10 m/s, which is to be 
compared with the experimental value of the sliding to 
stick-slip transition v c ~ 10~ 6 m/s. This very large dis- 



v b 

fs « (4-4.1)- 
and for Ni = 5, 
v bR ps 2.6-10- 2 . 



10- 3 , 

= 3, 
10- 2 ; 
while 



crepancy will be discussed in section [V. On the other 
hand, the friction fj, — |///i oa d| m the solid-sliding regime 
may be as small as fi ~ 10 -4 — 10~ 2 . 

Note that the value of v b in Eq. (|35| ) depends on the 
substrate mass ms which, in principle, may be taken ar- 
bitrary large. However, in fact only one or few boundary 
layers of the substrate play a role, and the results pre- 
sented above remains valid, at least qualitatively, for the 
infinite substrates as well P3j27|l . 



D. Friction-induced melting 

When the driving force approaches the threshold value 
/ = // , the energy pumped into the lubricant due to slid- 
ing can no longer be removed from the contact region into 
the substrates, and the lubricant temperature strongly 
rises up to its melting temperature. The solid-sliding 
regime is destroyed. To study this process, we started 
from the / = 0.2 steady-state configuration of the five- 
layer film, which is close to the threshold value // ps 0.31, 
and applied the dc force / = 0.5 to speed up the transi- 
tion. The results are shown in Fig. [L2[ The velocity v x 
of the top substrate grows approximately linearly with 
time. Because the top rigid substrate with one substrate 
layer attached to it moves as a whole during this process, 
we can use Newton's equation 

(m s + m S )v a = f^ ] = /M"--) + /(-ternal) ^ (gg) 

which for the x coordinate takes the form 2v x = 0.5 — /fries 
where /f r ; c is the frictional force acting on the top sub- 
strate from the lubricant. One can see from Fig. p"2| b 
that /f r i c ps 0.3 is the maximum force that lubricant 
losses can accommodate. Similarly one may find the 
normal force fp acting on the top substrate from the 
lubricant, 2v z = /i oa d + fp, which yields fp ps 0.1 (so 
that /i oa( j + fp ps 0) at the beginning, but at t > 20 
the force fp begins to grow and the width of the lubri- 
cant begins to increase. The lubricant temperature (see 
Fig. He) which was equal to T luh ps 0.2 for the / = 0.2 
steady state, at t > 20 begins to grow up to Ti UD ~ 0.6 
- 0.8 which is much higher than the melting tempera- 
ture T me it ~ 0.4, so the lubricant melts. Details of the 
forced transition are shown in Fig. |l^, where we plot the 
z coordinates of all atoms in the system. For t < 20 the 
lubricant slides as an ideal solid as it was for the / = 0.2 
steady state. At t ~ 20 the upper lubricant layers be- 
gin to slide over the lower layers, although the layers still 
keep the ideal triangular structure. At t ~ 30 the layers 
start to mix. At t ~ 40 the lubricant is already in the 
melted state, and the distribution of velocities across the 
lubricant v(z) exhibits a gradient. After the melting, the 
lubricant cannot anymore support the applied dc force 
/ = 0.5, and v top —> oo. 
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FIG. 12. "Forced" transition: the dependences of (a) the 
x velocities of the top substrate and the lubricant, (b) the 
normal and frictional forces, and (c) the lubricant tempera- 
tures versus time, when the force / = 0.5 is applied to an 
initial configuration corresponded to the / = 0.2 steady state 
of the five-layer film at T = 0. 
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FIG. 13. The z coordinates of all atoms and the top rigid 
substrate during the "forced" transition shown in Fig. [H| 



In another simulation we turned off the force at time 
t = 10.25, when the lubricant was still in the solid-sliding 
state, and allowed the system to relax. This situation 
mimics that of real systems, where the driving force exists 
due to elastic stress within the substrates. Thus, when 
the relative velocities of the substrate and the lubricant 
sharply increase, the stress is relaxed and the driving 
force decreases (in principle this situation can easily be 
modeled with a spring attached to the top substrate [Q; 
unfortunately, one has to introduce simultaneously few 
more model parameters) . In simulation we observed that 
although / = at t > 10.25, the top substrate continues 
to slide owing to its inertia. During this sliding the lubri- 
cant again melts as one can see from Fig. |l4[ Then the 
lubricant temperature decreases, and the lubricant solid- 
ifies into a six-layer (amorphous) film at t > 300 when 
Tiub ~ 0-3, and finally stops (note that the final configu- 
ration obtained with this force-induced melting and then 
freezing, is similar to the configuration obtained by in- 
creasing temperature till melting and then sharp decreas- 
ing it back to zero without the driving). The frictional 
force is large at the beginning when the lubricant melts 
(/fric ~ 0.15), but soon it quickly decreases to a value 

/fric ~ 0.01. 
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FIG. 14. The z coordinates of all atoms, when the driving 
force is removed at the time moment t = 10.25 of the "forced" 
transition of Fig. |l3| . 

We also made similar experiments but with an ap- 
plied force / ^ (/ < 0.5) at time t = 10.25. We ob- 
served that, if / < 0.005, the system stops in a six-layer 
state. For 0.005 < / < 0.03 we observed a sliding of the 
"frozen" six-layer lubricant; the lubricant temperature 
during the sliding is T lub w 0.06 at / = 0.005, T lub ~ 0.14 
at / = 0.01, and T lub « 0.19 at / = 0.02. For / = 0.03 
we observed that v — > oo, and the lubricant severs in the 
middle. From the slope of the v t0 p{t) dependences we 
found that /f r i c ~ 0.028, thus the frozen lubricant cannot 
support a driving force larger than / w 0.03. 

Comparing the results obtained for a realistic external 
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damping with the ones described above in the introduc- 
tion, where the model with a constant damping coeffi- 
cient was used, we see the essential difference. First, 
in the latter case the lubricant melts immediately when 
it begins to slide, i.e. ff < f s . This difference is due 
to the parameters chosen for the interatomic interaction. 
When the interaction of the lubricant with the substrates 
is much larger than the interaction within the lubricant, 
the scenario reviewed in the introduction is observed (the 
same also occurs if the lubricant is commensurate and 
perfectly aligned with the substrates). Second, after the 
shear-induced melting, the layer-over-layer sliding was 
observed, while now the transition to a "gas" phase takes 
place. Third, the back transition with force decrease now 
proceeds into an amorphous (frozen) state and not into 
the ideal solid one. These last two differences are due to 
damping mechanism used in the simulation. As will be 
shown below in Sec. IV, when the model with a constant 



damping is used, the scenarios of the previous works are 
observed, which points out the importance of the damp- 
ing mechanism on some results. 



E. Sliding of curved substrates 

The results presented above, were obtained for ide- 
ally flat substrates. Due to periodic boundary condi- 
tions we have in fact infinite substrates and the lubricant 
atoms are strictly confined between them. In a real sys- 
tem, however, the lubricant may leave the contact region 
through open boundaries. Besides, the substrates may be 
not ideally flat. To study this situation, we made a few 
simulations for curved substrates for the Ni = 5 system. 
The results are the following. 

The v top (/) dependence for the case when only one of 
the substrates (e.g., the top one) is curved in the x direc- 
tion (along the driving direction) is presented in Fig. [l5| 
by up triangles. In this case the lubricant film is approx- 
imately flat and follows the bottom (flat) substrate. At 
very low forces, / < 5d0 -4 , the lubricant moves together 
with the top substrate. For forces within the interval 
5 • 10~ 4 < / < 0.01 the ideal solid-sliding regime ex- 
ists. The lubricant film is approximately flat and slides 
over the flat bottom substrate, and the top substrate 
slides over the lubricant; in the contact region the lubri- 
cant is slightly compressed. For forces within the interval 
10~ 3 < / < 10~ 2 the velocity is approximately two times 
smaller than it was for the flat substrates. Note that be- 
cause the real contact area takes now only half of the 
total surface, the real load pressure in the contact is two 
times larger. At higher forces, / > 10 -2 , the perturba- 
tion of the lubricant becomes significant. At / > 0.011 
some atoms from the topmost lubricant layer begin to 
escape into the empty space between the lubricant film 
and the curved top substrate, although they are pushed 
back when the bottom of the curved top substrate reaches 
them during the sliding. This mixing increases with in- 



creasing of the dc force, and at / > 0.02 the lubricant 
fills the empty space. In the thick part there are six lay- 
ers. Finally at / sw 0.05 (i.e. much earlier than for the 
flat substrates) the lubricant melts, and the solid-sliding 
rea:ime is destroved. 
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FIG. 15. The x velocity of the top substrate versus force 
for the five-layer lubricant film when the top rigid substrate 
is curved (^ up) = 1, h^ p) = hi dn) = h[ dn) = 0, the up 
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When both substrates are curved in the x direction 
with the same curvature parameters, the ground state 
of the system corresponds to the five-layer lubricant film 
with an ideal structure of the layers. The curvatures of 
the substrates coincide, and the lubricant film just follows 
the substrates. But when the top substrate moves, one 
have alternatively a configuration where the curvatures 
coincide as in the ground state, and a saddle configura- 
tion where the curvatures anti-coincide so that in one half 
of the contact area the film is compressed, while in the 
other half the distance between the substrates is much 
larger than the width of the five-layer lubricant film. Be- 
cause a thin lubricant film cannot stay free in space but is 
attracted to one of the substrates, it has to follow cither 
the bottom or top substrate. In simulation we observed 
that the lubricant indeed takes a ground state-like con- 
figuration when the substrate curvatures coincide, while 
in the saddle points the lubricant follows the top sub- 
strate at the left-hand side, where the z-coordinate of 
the top substrate goes up (so it "drags" the lubricant), 
and at the right-hand side, where the z-coordinate of the 
bottom substrate goes up, the lubricant follows the bot- 
tom substrate. Thus, in the middle point the lubricant 
has to switch between these two configurations, and this 
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decreases the system mobility. The Vt op (f) dependence 
for this case is shown in Fig. [l5| by down triangles. The 
backward transition now takes place at /& « (2 — 3)-10 -3 . 
When the force increases up to / ~ 0.04, the film does 
not have enough time to reach the five-layer ground state- 
like configuration between the saddle points, ft first takes 
a six-layer structure but soon (at // > 0.04) the lubri- 
cant melts and the sliding is destroyed. Note, however, 
that the situation described above corresponds to a rare 
case when the curvature parameters of both substrates 
are exactly equal one another. In all other cases, in a 
given region only one of the substrates is usually strongly 
curved. 

We also checked the system mobility when the sub- 
strates are curved in the y direction (perpendicular to 
the sliding direction). When only the top substrate is 
curved, the velocity remains approximately the same as 
that for the flat substrates (for example, for / = 2-10 -3 
the velocity decreases from v « 0.98 for the flat sub- 
strates to Wtop ~ 0.96 for the curved top substrate). 
When both substrates are curved, the velocity remains 
unchanged (e.g., for / = 2-10~ 3 the velocity increases to 
"ftop ~ 0.985). Finally, when the top substrate is curved 
in both x and y directions simultaneously, the velocity de- 
creases for the / = 2-10 -3 case to the value vt op ~ 0.52 
(compare with v top ~ 0.57 for the case when the top 
substrate is curved in x direction only). 

Thus, we may conclude that the simulation results ob- 
tained for the system with flat substrates stay valid, at 
least qualitatively, for curved ones. 



IV. DISCUSSION AND CONCLUSIONS 

To summarize, we have studied the friction for a thin 
solid lubricant film between two flat solid substrates and 
found that the frictional losses are mainly due to the exci- 
tation of z vibrations in the lubricant with the washboard 
frequency. The anharmonicity of lubricant vibrations is 
found to be important and leads to heating of the lubri- 
cant. Thus, a first lesson learned from the simulation of 
a complex model, is that simple FK-type models with 
harmonic springs between the atoms cannot describe the 
kinetic friction properly. 

The solid lubricant between two flat (or at least not too 
rough) substrates provides the smallest possible frictional 
force for high-speed systems. The energy losses are very 
small both in the low-velocity (but Vt Qp > Vb ~ 10 m/s) 
case, as well as in the high-velocity (vt op ^ fsound) 
regime, the later regime being stable if there is a gap 
between the lubricant and substrate phonon spectra. If 
the lubricant film is in a liquid state (melted either due 
to high substrate temperature or because of very high 
speed), the frictional force is much larger as seen from 
Fig. O. When the lubricant is frozen back from the 
melted state, it takes an amorphous structure because, 
due to the contact with the substrates, the energy is re- 



moved from the lubricant very fast, and the confined lu- 
bricant film has no time to order. The frictional force in 
the case of an amorphous lubricant may be even larger 
than for the melted one (see Fig. rLq). 
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FIG. 16. The v t0 p(f) dependences for the melted lubricant 
(Tiub ~ 0.5, the up triangles) and the frozen (amorphous) 
lubricant (the down triangles). 

The solid lubricant between flat substrates leads also 
to quite small values of the static frictional force f s . This 
agrees with the results of previous simulations (see and 
references therein). The new result of the present work 
is the dependence of f s on the number of layers in the 
solid lubricant film. The exponential decrease of f s with 
Ni may be used for explanation of experimental values 
of the velocity of the transition from smooth sliding to 
stick-slip motion as will be discussed below. 

We also checked how the results are sensitive to the 
parameters of the model. First, changing the mass ms 
of the top rigid substrate practically does not change 
the adiabatic Vt op (f) dependence, thus the velocity in 
the steady solid-sliding regime does not depend on ms- 
Another important parameter is the minimal external 
damping 77 m j n which models the multi-phonon and/or 
"electron-hole" contributions to the energy exchange of 
lubricant with substrates. We found that the results are 
approximately unchanged for vt op 1, when the energy 
losses inside the lubricant play the main role. At small 
forces/velocities, when vt op ^ 1, the values of Vt op and 
Zto P are larger for smaller values of ?7 m i n and agree with 
the "perfect-sliding" dependence (|25|). 

Finally, let us discuss the relation of the microscopic 
simulations of the present work to macroscopic friction. 

Dependence on load. We varied /i oa d to test the Amon- 
tons first law Ff r i C = /i-Fioad, where y, (called the friction 
coefficient by tribologists) should be approximately con- 
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stant. For example, the simulation results for a much 
lower load /i oa d = — 10~ 4 than that used in all simula- 
tion presented above, are shown in Fig. |l7|. 
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FIG. 17. The dependences Vt op (f) for two values of the 
load /ioad = —0.1 (dashed curve and solid diamonds) and 
/load = — 10~ 4 (solid curve and open diamonds) for the 
five-layer film at T = 0. Inset: «to P (/ioad) (left axes) and 
2to P (/ioad) (right axes) for the fixed value of the driving force 
/ = 0.02. 



ulations. Besides, the authors of [|19| used the Langevin 
equations with a constant external damping. 

Dependence of results on the interaction parameters. 
Above we have studied the system behavior for the fixed 
set of parameters, when the interaction within the lubri- 
cant, Vu = 1, is much larger than the interaction of the 
lubricant with substrates, V s i — 1/3. For these parame- 
ters the solid-sliding regime had indeed to be expected. 
Below we present a few simulation results for other sets 
of narameters. 
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The trivial result is that the width of lubricant in- 
creases with decreasing of /i oa d- Therefore, the external 
damping coefficient, which exponentially depends on z, 
decreases with /i oa d- As a result the backward transition 
now takes place at a lower force / = /b ~ (3 — 4) ■ 1CP 5 
when Vb < 3-10~ 2 . The forward transition also occurs 
earlier, at ff ps 0.1 — 0.11 when Vb 3. For v t0 p ^ 1 t ne 
dependence Vt op (f) is approximately the same as that for 
the high load, while for lower forces when v top < 1, the 
values of vt op for the same values of / are larger than 
those for the /i oa d = —0.1 case, because the external fric- 
tion is smaller due to larger values of z. The dependences 
of «top and Ztop on /i oa d for a fixed value of the dc force 
/ = 0.02 are shown in inset of Fig. [l7|. Surprisingly, ti top 
is not decreasing but increases when /i oa d increases, be- 
cause for the chosen value of the dc force, a resonance 
effect plays a role (the washboard frequency penetrates 
into the lubricant phonon spectrum). In general, how- 
ever, the variation of /i oa d results in a little change only. 
This confirms once more that the empirical Amontons 
first law works due to change of the contact area with 
load H. On the contrary, He and Robbins jl9| have 
found that the frictional force is directly proportional to 
the load. However, their simulation has been done with 
the const-D algorithm at low velocities (ftop < 0.1 n.u.) 
and much higher temperatures than studied in our sim- 



FIG. 18. The «top(/) dependence for the one-layer lubri- 
cant with the interatomic interactions Vu = 1/3 (solid curve 
and open diamonds) comparing with that for the Vu = 1 case 
(dashed curve and solid diamonds). Inset: the Vu = 1/3 (solid 
curves) and Vu = 1 (dotted curves) spectra for the five-layer 
lubricant. 

First, let us consider the case with Vu = V s i = 1/3, 
when the interaction within the lubricant has the same 
magnitude as the lubricant-substrate interaction. In 
fact, this case is similar to that studied above, because an 
effective lubricant-substrate interaction is again smaller 
than the interaction within the lubricant due to mismatch 
of the substrate/lubricant lattice constants. However, an 
essential new feature of the present case is that now the 
phonon frequencies of the lubricant are smaller. Because 
the interaction within the substrates and within the lu- 
bricant differs by nine times, the substrate and lubri- 
cant phonon spectra are well separated one from another 
even for the five-layer lubricant film as shown, for exam- 
ple, in inset of Fig. |l8| for the z component. There is 
a wide gap for the frequencies 4 < uj < 6.5 where the 
density of phonon states is very small. Thus, one could 
expect that if the washboard frequency comes into this 
frequency interval, the effective friction should be small, 
and the solid sliding regime will persist up to very high 
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velocities. Indeed, the simulation results for the one-layer 
lubricant film presented in Fig. [l8], show that the solid- 
sliding regime indeed survives now till the very high force 
(// ~ 5) and velocity (i 1 / ~ 10). This effect is similar 
to the "reentrant" transition observed by Stevens and 
Robbins Jl6|], where with the increase of the driving ve- 
locity the lubricant melts but then, at a higher velocity, 
it solidifies again. In our simulation with a small ex- 
ternal damping we did not obtained melting but instead 
observed a large increase of the lubricant temperature. 
For velocities v top > 3 (recall that in the previous case 
the solid-sliding regime was already destroyed at so high 
velocities) the washboard frequency is higher than the lu- 
bricant phonon frequencies, thus the excitation of oscil- 
lations in the lubricant as well as the effective lubricant 
temperature decrease. This regime corresponds to the 
"flying" one, i.e. the decoupling of phonon spectra of the 
lubricant and the substrates stabilizes the "flying" state. 
With increasing dc force the lubricant remains sliding in 
a solid state, while the amplitude of the z oscillations 
of the substrate atoms strongly increases at u top > 6, 
and this finally destroys the solid-sliding regime. Note 
that for the forces 0.1 </ < 0.4, when 2 < v t0 p < 3, 
the system behavior is similar to that described above, 
because the densities of lubricant vibrations are simi- 
lar at these washboard frequencies. At smaller forces, 
/ < 0.1, the effective friction as well as the lubricant 
temperature are larger for the weak-lubricant case than 
those for the previous case. The backward transition to 
the locked state takes place now at ft < 0.003, when 
v b w 0.15. The static frictional force, f s w 0.013-0.014, 
is also smaller than that for the previous case (recall 
f s « 0.02 for the Vu — 1 case). Note, however, that 
we observed the "flying" regime only for the one-layer 
film. For the five-layer lubricant the forward transition 
takes place at // « 0.025 — 0.026, i.e. much earlier than 
for the Vu = 1 case, and the attempts to get a high- 
velocity regime failed. Note that although it seems quite 
unusual to have a velocity close to the sound speed in 
real machines, the "flying" regime may be reached in, 
e.g., next generation of hard-disk drives. In a disk with 
a radius of 4 cm and a rotating speed of 10000 rpm, the 
speed of the head with respect to the plateau already ex- 
ceeds 40 m/s, i.e. the crash of the head on the plateau in 
a failure is an example of high speed friction, although, 
unfortunately, the "flying" regime would not be reached 
the current generation of hard disks. 

Second, we studied the opposite case, when the 
lubricant-substrate interaction is much larger than the 
interaction within the lubricant, namely V s i = 3 and 
Vu = 1/3. In this case the two utmost lubricant lay- 
ers tend to be attached to the corresponding substrates, 
so a sliding can only emerge inside the lubricant. In fact, 
this situation is close to real oil lubricant, where the in- 
teraction with substrates prevents the lubricant from be- 
ing squeezed out from the contact areas. In the case of 
Ni = 2 all lubricant atoms are attached by the substrates, 
so there are no "lubricant" between the two sliding "bod- 



ies" , each consisting of the rigid substrate with one at- 
tached s-layer and one "glued" layer of 1-atoms. Thus, 
the situation is close to the dry friction for two commen- 
surate perfectly aligned substrates. However, the interac- 
tion between the bodies, Vu — 1/3, is much smaller than 
the interaction inside them, where V ss — V s i — 3. In this 
case at some critical force, f s ~ 0.08, the bodies decouple 
one from another, the distance between them increases, 
and they begin to "fly" almost without friction. But 
the loading force pushes the bodies into contact again, 
so they collide and almost elastically repulse, then again 
collide, etc. The motion of the top substrate resembles a 
series of jumps as shown in Fig. [1^. If the driving force 
remains constant, the amplitude of jumps increases with 
time, and finally at one collision the surface lubricant 
layers are destroyed. 



15 



10 



5 - 



Vdeft) 

z , op fag 111 ) 

V=l/3 

V sl =3 
V =3 

ss 

N=2 



/ 

A 

I A i 1 

/ i j 1 



/4 



v i 



A; w « i! 



80 



60 



top 



40 



20 



200 



400 
time 



600 



800 



FIG. 19. Beginning of motion for the "glued" two-layer 
lubricant film (V as = Vi s = 3, Vu = 1/3). 

For a thicker lubricant film, Ni > 2, there are weakly 
interacting lubricant atoms between the two bodies. Note 
that for these parameters in the ground state configura- 
tion the utmost lubricant layers which are attached to 
the substrates, have slightly more atoms than the middle 
layer (s). When the motion begins at / = f s in this case 
(e.g., for the N t = 3 or N t = 5 system f s ~ 0.03 - 0.07), 
the middle lubricant layers melt according to a 3D melt- 
ing mechanism (one layer in the Ni = 3 system or three 
layers for the Ni — 5 case) , and the width of the lubricant 
film increases. Starting with the configuration just after 
the beginning of the motion (e.g., that for / = 0.08 for 
the Ni = 5 system), we abruptly decreased the dc force 
to a smaller value. The results are presented in Fig. ||(]. 
If / < 0.01, the motion locks in the five-layer solid state, 
while for / > 0.02 the top body is "flying" and its veloc- 
ity increases linearly with time. Thus, the frictional force 
for a melted lubricant in this system is /f r i c ~ 0.015. 
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FIG. 20. Motion of the "glued" five-layer lubricant film 
(Vss = Vis — 3, Vu = 1/3) for different driving forces. 

Third, all the scenarios described above were obtained 
for the model where the external damping depends on 
the distance from the substrates and, thus, it is very 
small for the middle lubricant layer (s). However, if we 
change the model in the way that the minimal external 
damping would not depend neither on distance nor on 
velocity, the scenario changes to that already observed 
in the previous works []l5|,[nj and described in the In- 
troduction. Now when the sliding begins the width does 
not increase significantly and the sliding steady state sur- 
vives with increasing the dc force. With increasing force, 
the middle layers melt due to sliding according to a 2D 
melting mechanism, i.e. the lubricant keeps its layered 
structure and different lubricant layers slide one over an- 
other. This model may be applied to lubricants with 
complex molecules like organic ones, where an internal 
viscosity of the lubricant plays an important role. 

The transition from smooth sliding to stick-slip mo- 
tion. If, instead of a constant dc force we would drive 
the system through a spring attached to a stage mov- 
ing with the velocity vtop, the system would exhibit a 
transition from smooth sliding to stick-slip motion at 
the velocity «top ~ Vb f27|]. A similar result would be 
expected with a large enough substrate block so that 
its elasticity would start playing a role. Experimentally 
the transition from smooth sliding to stick-slip motion 
is observed in almost all systems where the static fric- 
tion force is nonzero (the zero f s corresponds to liquid 
lubricant). There are at least three mechanisms that can 
explain the transition from smooth to stick-slip motion: 
(a) "inertia" mechanism, when the backward sliding-to- 
locked transition occurs at the force ft < f s due to in- 
ertia of the moving lubricant in the underdamped sys- 



tem (this mechanism was discussed above in Sec. Ill C ) ; 
(b) "melting-freezing" mechanism |]l5|,p|,p|jic|l , when the 
lubricant undergo dynamical phase transitions between 
a fiuidized state during slip (sliding) and a solid state 
during stick; and (c) "memory" mechanism (e.g., see 
l@,^|Jl])i when, after the sliding-to- locked transition, the 
static friction force increases with the time of station- 
ary contact. The simulation results of the present work 
show that the "inertia" mechanism leads to the criti- 
cal velocity v c ~ 10~ 2 n.u. w 10 m/s, and we do not 
see any reason to get a much lower value for v c . This 
value is orders of magnitude larger than the experimen- 
tal value, which rules out the inertia mechanism to ex- 
plain the experimentally observed stick-slip to smooth 
sliding transition. For the model parameters used in the 
present work, we also did not observed the melting of 
thin lubricant films except for extremely large forces and 
velocities. The lubricant may melt at / > f s if it is 
strongly coupled with the substrates or, in a more spe- 
cific case, when the lubricant is perfectly commensurate 
and aligned with the substrates |l5|]27| ]. But, in such sit- 
uations, when the driving force is reduced, the lubricant 
quickly solidifies again due to a fast energy exchange with 
the substrates (as discussed above in Sec. Ill D| ) , and the 
backward transition to the locked state takes place at the 
velocity v c ~ 1 n.u. [ ^5||3^j27[] . Thus, in this case too the 
top substrate cannot move smoothly with a low velocity 
at the microscopic scale. In experiments the transition 
from stick-slip to smooth motion is observed usually at 
a velocity v c ~ 1 fxm/s ~ 10~ 9 n.u. [Q. Even if we 
suppose that there are only few contacts (asperities) be- 
tween the solids and only one of the contacts moves at 
a time, we do not see a way by which the macroscopic 
stick-slip to smooth transition could be explained by the 
microscopic (atomic-scale) mechanisms. Thus, we con- 
clude that the mechanism responsible for the experimen- 
tally observed stick-slip to smooth transition, at least for 
"simple" lubricant molecules, should be the "memory" 
effects. After the sliding-to-locked transition, the static 
friction force increases with time due to plastic deforma- 
tion of the lubricant, for example, because of an increase 
of the area of real contact, or due to the decrease of the 
lubricant width (according to the simulation results of 
the present work, the squeezing of the lubricant leads to 
the increase of / s ), or due to interdiffusion of lubricant 
molecules between different layers of the lubricant in the 
case of lubricant consisting of long-chain polymers. All 
these processes involve plastic deformations and thus are 
characterised by macroscopic-scale characteristic times. 
Then, the experimentally observed velocities of the tran- 
sition from smooth sliding to stick-stick motion can be ex- 
plained with the help of earthquakes-like models ]3r^-|3^1 
as was demonstrated in p j39| , fIo|,p7| | . Consequently the 
analysis of the order of magnitude of the stick-slip to 
smooth-sliding transition velocity suggests that the ex- 
planation of this penomcnon should not be seeked at the 
microscopic scale and is probably out of reach of current 
molecular dynamics simulations. Besides, this also ex- 
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plains why the experimentally measured kinetic friction 
in the smooth sliding regime is often almost independent 
of the velocity, while the microscopic one shows a very 
strong dependence. Indeed, if the macroscopic smooth 
sliding corresponds to the microscopic stick-slip regime, 
the macroscopic kinetic friction would correspond to mi- 
croscopic static friction force which is constant. There- 
fore the evaluation of the orders of magnitude that we 
made in this work lead us to conclude that there is a 
fundamental difference between microscopic friction and 
the macroscopic experimental properties. This should 
perhaps be considered for the current projects to build 
micro- machines. 

Finally, in real experimental systems the rate of change 
of the driving force corresponds to adiabatically slow 
variation of / in simulations. Therefore, to obtain reli- 
able results, which do not depend on an initial configura- 
tion, it is important to reach the steady state for a given 
/. As observed in the present work, lubricant systems 
are often characterized by quite long relaxation times 
Troiax > 10 3 n.u. « 10~ 10 s (or even T rc i ax ^> 1CP 10 s in 
the solid-sliding regime at velocities Vtop < 10 m/s). In 
simulations for large systems with realistic interactions 
one could expect to reach the stationary state for one 
given set of parameters (including the value of the driv- 
ing force) only, and not for adiabatically changing force. 
Thus, although being very attractive, the "millions of 
atoms" simulations hardly could give reliable values for 
the kinetic friction in their present stage. This is why 
we believe "simple" models are still very important. But 
they must also be taken with great caution because, as 
we have shown in the present work, the results obtained 
with simple models may also strongly depend on the 
model chosen, in particular of the damping introduced 
in the simulations. We tried here to introduce a damp- 
ing mechanism motivated by microscopic considerations 
and we showed that the results are sometimes very dif- 
ferent from those given by a uniform damping coefficient 
as often considered. 
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